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Abstract 

Anomalous behaviour of the excess entropy {Se) and the associated scahng relationship 
with difFusivity are compared in liquids with very different underlying interactions but similar 
water- like anomalies: water (SPC/E and TIP3P models), tetrahedral ionic melts (Si02 and 
BeF2) and a fluid with core-softened, two-scale ramp (2SRP) interactions. We demonstrate 
the presence of an excess entropy anomaly in the two water models. Using length and 
energy scales appropriate for onset of anomalous behaviour, the density range of the excess 
entropy anomaly is shown to be much narrower in water than in ionic melts or the 2SRP 
fluid. While the reduced diffusivities {D*) conform to the excess entropy scaling relation, 
D* = Aexp{aSe) for aU the systems (Y. Rosenfeld, Phys. Rev. A 1977, 15, 2545), the 
exponential scaling parameter, a, shows a small isochore-dependence in the case of water. 
Replacing Se by pair correlation-based approximants accentuates the isochore-dependence of 
the difFusivity scaling. Isochores with similar diffusivity scaling parameters are shown to have 
the temperature dependence of the corresponding entropic contribution. The relationship 
between diffusivity, excess entropy and pair correlation approximants to the excess entropy 
are very similar in all the tetrahedral liquids. 
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1 Introduction 

Water displays a number of thermodynamic and kinetic anomalies when compared to simple 
liquids [TH3]. The density anomaly, corresponding to a negative isobaric thermal expansion 
coefficient (ap), is the best known of these unusual properties of water and is observed for state 
points lying within an approximately parabolic boundary defined by the locus of temperatures 
of maximum density (TMD) for which ap = 0. The density anomaly implies the presence of 
other thermodynamic anomalies, such as those associated with the isobaric heat capacity (Cp) 
and the isothermal compressibility (kt)- The kinetic anomalies of water are associated with 
an increase in molecular mobility on isothermal compression, measured in terms of diffusivity, 
orientational relaxation times or viscosity. A number of network-forming inorganic melts with 
local tetrahedral order have been shown possess water-like anomalies, most notably, AB2 ionic 
melts such as Si02 and BeF2 and the liquid phase of elements, such as silicon and tellurium 
[3Vtl5). More recently, water-like anomalies have been demonstrated in mesoscopic liquids with 
isotropic, core-softened effective interactions or anisotropic patchy interactions |16fll8| . Despite 
a very diverse set of underlying interactions, liquids with water-like anomalies are found to 
have essentially conformal liquid-state "phase diagrams" with a nested structure of anomalous 
regimes of density, diffusivity and structural order. 

The similarity in the phase diagrams of water-like liquids reflects similar structure-entropy- 
diffusivity relationships that can be conveniently analysed in terms of the excess entropy, Se, 
defined as the difference between the total thermodynamic entropy (S) and the correspond- 
ing ideal gas entropy (Sid) at the same temperature and density. A necessary condition for 
a fluid to show water-like thermodynamic and transport anomalies is the existence of an ex- 
cess entropy anomaly, corresponding to a rise in excess entropy, Se, on isothermal compression 
{{dSe/ dp)T > 0) [9 HTH[T9H25] . Liquids with water-like anomalies display distinct forms of local 
order or length scales in the low- and high-density regimes; competition between the two types 
of local order results in a rise in excess entropy at intermediate densities. Most liquids, includ- 
ing anomalous ones, obey semi-quantitative excess entropy scaling relationships for transport 
properties of the form X* = Aexp{aSe) where X* are reduced transport coefficients, and A and 
a are scaling parameters that are very similar for systems with conformal potentials [26H29j . 
Consequently, the existence of an excess entropy anomaly is reflected in mobility anomalies. In 
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order to make a more precise connection between thermodynamic and mobility anomalies, it is 
necessary to understand the assumptions underlying the scaling rule. Rosenfeld-scaling assumes 
that diffusion in liquids takes place through a combination of binary collisions and cage relax- 
ation. The binary collision contribution is approximately factored out by using macroscopic 
reduction parameters based on elementary kinetic theory. For example, reduced diffusivities 
are defined as D* = D{p^/^ / (mkBT)^/'^) [27\. The frequency of cage relaxations is assumed to 
be proportional to the number of accessible configurations, exp(aS'e), since configuration space 
connectivity is high in the stable liquid phase. In order for the scaling relationship to be state- 
point independent, it is necessary that the exponential parameter a controlling the number of 
accessible configurations is determined by the interaction potential and is otherwise state-point 
independent. 

This paper develops a basis for quantitative comparison of structure-entropy-diffusivity rela- 
tionships in liquids with very different underlying interactions but similar water-like anomalies. 
We focus on three different categories of liquids: (i) molecular fluids (H2O), (ii) tetrahedral 
ionic melts (BeF2 and Si02) and a core-softened fluid with isotropic, pair-additive interactions. 
Comparison of these different systems should provide very useful physical insights into the re- 
lationship between structure, entropy and transport properties of liquids, including complex 
fluids, and has not been previously attempted. It is evident that parametric variations of a sin- 
gle functional form for the interaction potential are insufficient to serve as a basis for comparison 
of anomalous behaviour for the three different categories of liquids studied here. Instead, we 
focus on directly comparing the two crucial features of the excess entropy relevant for water-like 
liquids: firstly, the presence of an excess entropy anomaly and secondly, exponential scaling of 
transport properties with excess entropy. Since the state point, {pm,Tm), corresponding to the 
maximum temperature along the TMD locus is a unique point common to all liquids with a 
water-like density anomaly, we use this information to define a natural energy and length scale 
for onset of anomalous behaviour in each system. 

The specific interaction models used for the different categories of liquids require some 
introductory comments. Anomalous behaviour and scaling of transport properties with the 
thermodynamic excess entropy has so far not been studied in the case of water. Since the 
different water models are known to display anomalous properties at very different temperatures, 
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we use two different effective pair potentials for water, TIP3P and SPC/E. The SPC/E model is 
known to qualitatively reproduce all anomalous properties of water but typically at temperatures 
that are 30K to 40K lower than the experimental value of 279K [SHU [191 [30]. The TIP3P 
model is widely used in biomolecular simulations, but the onset temperature for anomalous 
behaviour is 170-180K |19 p 3 H [32]. The transferable rigid-ion model (TRIM) potential was used 
for BeF2 [33] and and the van Beest-Kramer-van Santen (BKS) potential was used for Si02 [34j . 
As an example of a core-softened fluid with waterlike anomalies, we consider the two-scale ramp 
potential (2SRP) [16pi7j . Thus we study five different simulation models for the three categories 
of liquids. Table 1 summarises some of the relevant information about these simulation systems 
at {pm,Tm)- Section 2 provides the necessary simulation details. Results are presented in 
Section 3 and conclusions are discussed in section 4. 

Simulation Details 

Potential energy models for all systems studied here are given in Table 1. Molecular dy- 
namics simulations of ionic melts (BeF2 and Si02) and water models (TIP3P and SPC/E) were 
performed in the NVT ensemble using the Verlet algorithm as implemented in the DL_POLY 
package [35]. Temperature range of 1500-3000K is studied for BeF2 and 4000-6000K for Si02. 
Details of simulations of the two ionic melts are given in ref. [36]. Simulations of TIP3P and 
SPC/E water used 256 molecules in a cubic simulation box with 1 fs timestep and production 
run lengths of 4 to 8 ns. Rigid-body constraints were maintained using the SHAKE algorithm. 
The Berendsen thermostat time constant, tb, is Ips for SPC/E for all state points except those 
along the 210K isotherm (rB=200ps) while tb= 20ps for all TIP3P state points. State points 
for SPC/E water covered a temperature range from 210K to 300K and density range from 0.9 
g cm~^ to 1.4 g cm~^; the temperature range for TIP3P water was kept as 170K to 300K. 

Simulations of the two-scale ramp potential (2SRP) fluid were performed using Metropolis 
Monte Carlo simulations using 256 particles in a cubic simulation box. The 2SRP fluid has 
length scales associated with the hard-core (do) and soft-core (fii) diameters; we use ai as the 
unit of length. The energy Ui of the linear ramp interaction extrapolated to zero separation is 
taken as the unit of energy. State points cover a temperature range of 0.025 to 0.2 and density 
range from 0.8 to 2.6 in reduced units j^). Diffusivities for five isotherms from T = 0.027 to 
T = 0.063 were taken from ref. j^37j. 



Widom insertion method was used to calculate chemical potential and excess entropy for the 
2-scale ramp fluid at low densities and high temperatures [38]. Thermodynamic integration was 
then used to obtain excess entropies at other statepoints. For ionic melts, Sid is computed for 
a non-interacting multi-component mixture of atoms with the same masses and mole fractions 
as the ionic melt. Excess entropy for BeF2 and Si02 melts were calculated using the same 
procedure as in ref [19], but over a more extensive density range. The excess entropy of the 
water models was computed relative to an ideal gas of rigid triatomic molecules at lOOOK, 
O.Olg/cc. Thermodynamic integration was used to compute Se at other state points [4J. Our 
entropies at 298K and 1 g cm~'^ match those of [39J for SPC/E and TIP3P, within statistical 
error. The use of atom-atom RDFs in equation(l) would suggest that an ideal gas reference state 
analogous to that used for ionic melts would be more appropriate; this is in fact not necessary, 
since the difference between the two reference states will be constant at a given temperature 
and density. We note that while the system sizes chosen were chosen to be small enough that a 
large number of state points could be efficiently covered, they are sufficiently big that finite-size 
effects would be expected to make only small quantitative differences to the computed values 
of thermodynamic and transport properties. 

3 Results 

3.1 Excess Entropy Anomaly 

The strength of the excess entropy anomaly determines the presence of diffusional, density 
and related anomalies and reflects the variation in structural order in a fluid as a function of 
density. Figure 1 compares the density dependence of the excess entropy along the T = Tm 
isotherm, scaling densities with respect to the pm value of each system. While the overall 
behaviour of Se{p/pm) is very similar, it is immediately apparent that, in units of /Om, the 
water models have a relatively narrow density regime of anomalous behaviour and a steeper 
slope in the high density regime. The two-scale ramp (2SRP) system in contrast, has a much 
wider anomalous regime bounded by a well-defined minimum as well as maximum in Se{p) 
along the Tm isotherm. The two ionic melts resemble the ramp fluid in the density range of 
anomalous behaviour and the location of the maximum in Se{p/pm)- The minimum in the 
Se{p/ Pm) curves could not, however, be located within the density range consistent with the 
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thermodynamic stability criterion that the isothermal compressibility should be positive. In 
units of Nks, the rise in the Se{p) curves from the minimum to the maximum is a fraction of a 
ks for all the systems. This small anomaly is, however, sufficient to create an inflexion point in 
the total S{p) = Sid + 'S'e curves that marks the onset of the density anomaly. The qualitative 
behaviour of the excess entropy anomaly is consistent with the grouping of the fluids into three 
similar classes: rigid water models, ionic melts and core-softened fluids. Unless otherwise stated, 
results for only one member of each class will be discussed in the rest of this paper. 

3.2 Excess Entropy Scaling of DifFusivities 

A comparison of the excess entropy scaling of the diffusivities for the different systems 
studied in this work, shown in Figure 2, illustrates that such a scaling approach is a reasonable 
predictor of diffusivities though there are small system-dependent differences. The 2SRP fluid 
shows excellent scaling for the state points studied, but the range of variation of the diffusivities 
is small in comparison to both ionic melts and water, given temperature variations by factors 
of two to three. Both Si02 and BeF2 show very good scaling in the liquid regime, but there is 
a change in slope a for the lowest temperature isotherms, possibly due to onset of cooperative 
dynamics |21il221 [36 l l40 p 41j . The water models all show a distinct isochore dependence of the 
In D* versus Se plots, which is less pronounced in ionic melts and is essentially zero in the 2SRP 
fluid. In the case of water, deviations from linearity are specially noticeable at low temperatures 
for the 0.90 and 1.40 g cm~^ isochores for which the mode-coupling temperature is relatively 
high [42j. All liquids with water-like anomalies are associated with a competition between two 
different length scales or local order metrics as a function of density. Therefore deviations from 
corresponding states behaviour would be expected. It appears, however, that these deviations 
are small and most noticeable in the case of water which has a rigid, non-spherical structure. 
Possible reasons for this are discussed below. 

3.3 Pair Correlation Estimators of the Excess Entropy 

The effect of the underlying structural correlations in a fluid on the entropy-transport 
relationships can be understood by expressing Se as a multiparticle correlation expansion, 
Se = S2 + S3 + . . ., where Sn is the entropy contribution due to n-particle spatial correla- 
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tions [43fl45] . The entropy contribution due to pair correlations between atoms of type a and 
P is given by: 

Sap= f {gai3{r)\ugap{r) -[gap{r) -l]}r'^dr (1) 
Jo 

where Qapif) is the atom-atom pair distribution function (PDF). The overall pair correlation 
entropy, 52, is: S2/NkB = -'^'^pY.a. 

^ XiyXj^jSafj where A'^ is the number of particles and Xq, 
is the mole fraction of component a in the mixture. 5*2 can be estimated from experimental 
scattering data as well as simulations and is the dominant contribution to the excess entropy 
in the case of simple [Hj and core-softened fluids (e.g. results for 2SRP fluid in this paper), 
ionic melts and water [19]. In the case of the AB2 melts and water, we also compute the 
tetrahedral correlation entropy, Saa, due to the pair correlations between the tetrahedral(A) 
atoms, corresponding to Be, Si and O in BeF2, Si02 and H2O respectively. The Sqo entropic 
contribution is reproduced by many coarse-grained model potentials of water that replace each 
water molecule by a structureless particle with isotropic effective interactions that reproduce 
the behaviour of goo{f) (36t lT7]. 

3.3.1 DifTusivity Scaling 

The effect of using a structure-based approximant to the excess entropy, such as 52, on the 
Rosenfeld-scaling relation can be seen by substituting 5e = 52 + A5 to obtain 

D* ^ Aexp(a5e) = ^exp [a52(l + A5/52)] = ^exp(a'52) (2) 

where the new exponential scaling parameter, a' = I+AS/S2, has a state-point dependence that 
is determined by the density and temperature-dependent behaviour of A5/52. Thus by using 
different approximants to the excess entropy and observing the effect on the scaling behaviour 
of the transport properties, one can determine importance of different structural correlations 
in controlling mobility. Figure 3 shows the correlation between the dimensionless diffusivity 
(D*) and 52 as well as XqqSqo for SPC/E water, which should be compared with Figure 
2(a). It is immediately obvious that the temperature dependence along an isochore is linear 
or quasi-linear, but there is a distinct isochore dependence of the slope a of the InD* versus 
5^ plots. The isochore-dependence is most pronounced for 5^ = Sqo and least pronounced for 
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= Se- While the InD* versus S2 or Saa behaviour has been studied before [TSIESIESIES] , 
no systematic comparison with Se and 5*2 scahng has been made so far. The ionic melts show 
a parallel behaviour to the water models, in terms of diffusivity scaling with respect to Sg, 
S2 and Saa while the 2SRP fluid also shows a small isochore dependence when D* is scaled 
with respect to 82- In the case of BeF2, the diffusivity and viscosity scaling with respect to Se 
and ^2 is given in refs. [12l[36] and with respect to SseBe in Supporting Information Figures 
1. Scaling with respect to both Se and S2 for 2SRP fluid is shown in Supporting Information 
Figure 2. This is in contrast to simple liquids with pair additive potentials for where excess 
entropy scaling relationships based on Se and 52 have a similar degree of universality 128] , 

3.3.2 Temperature Dependence 

The isochore-dependence of diffusivity scaling with different entropic contributions, and the 
contrast with simple liquids, suggests that the isochoric temperature dependence of different 
entropic measures determines the quality of entropy-diffusivity scaling. Based on measure-free 
density functional theory for repulsive potentials, Rosenfeld predicted a T~'^l^ scaling for the 
excess entropy [H]. In the case of simple liquids, these predictions are satisfied by both Se and 
5*2 and Rosenfeld-scaling holds very well [481149] . Se is known to show scaling in ionic 

melts and water over temperature regime of interest in this study [5','6]. Figure 4 shows Se, 
S2 and 5*00 as a function of T"^/^ for different isochores in SPC/E water. Unlike in simple 
liquids, S2 and Sqo show significant deviations from T~'^l^ scaling in the anomalous regime. 
The results for BeF2 are very similar to those of SPC/E water (see Supporting Information 
Figure 3 and ref. [I0l[19]). In the case of the 2SRP fluid, both Se and S2 show deviations from 
7"-2/5 ijgiiaviour at low temperatures (see Supporting Information Figure 4). 

A comparison of the isochore-dependence of the Rosenfeld-scaling parameter with respect to 
an entropic contribution, 5^, and the temperature dependence of 5^ shows that isochores with 
very different temperature dependences also have very different diffusivity scaling parameters. 
For example. Figure 3(b) shows when InZ?* is plotted as a function of Sqo^ then the slope 
increases with increasing density and the three densest isochores, at /? = 1.2, 1.3 and 1.4 g cm~^ 
show very different and progressively increasing values of a. Figure 4(c) shows that Sqo along 
these three isochores has a distinctly different and very weak temperature dependence. The two 
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AB2 ionic melts show a very similar pattern of behaviour (see Supporting Information Figures 
1 and 3), reflecting the similar structure of the tetrahedral fluids. The 2SRP fluid shows a weak 
isochore dependence of a when is replaced by 5*2, which may be related to a wider spread of 
^2 values compared with S^- It is interesting to note that S2{T) and SaaIT) curves show a much 
stronger density dependence than Se(T). In the case of the tetrahedral liquids, it is notable 
that at high densities, Saa remains almost constant with temperature while diffusivity rises 
quite sharply. Since the overall excess entropy scaling for high and low density isotherms is not 
qualitatively different, this must imply that significant reorganisation of the network associated 
with the other pair correlation contributions contributes to controlling diffusivity; for example, 
in water this could be strong librational-translational coupling reflected in the Sqh and Shh 
contributions. 

4 Discussion and Conclusions 

The excess entropy anomaly is compared in five liquids with water-like anomalies, with 
three completely different types of interactions: ionic melts (BeF2 and Si02), hydrogen-bonded 
molecular liquid (SPC/E and TIP3P water), core-softened fluid (two-scale ramp). Using the 
state point, {pm,Tm), corresponding to the maximum temperature along the TMD locus to 
define energy and length scales, we show that excess entropy anomaly in the ionic melts and 
the core-softened fluid is very similar in range and strength, while the two water models have a 
conspicuously smaller density range and a much sharper decrease in Sg with density on strong 
compression. 

We show that scaling of diffusivities with the thermodynamic excess entropy is good in 
all cases, noting that water is the first molecular fluid, other than the Lennard -Jones chain 
fluid [50] , for which Rosenfeld-scaling with the thermodynamic excess entropy, rather than pair- 
correlation based estimators, has been tested. In comparison to ionic melts and the 2SRP fluid, 
a small isochoric dependence of the Rosenfeld-scaling parameters is seen in the case of water. 
This suggests the presence of additional length scales due to the rigid, non-spherical molecular 
shape, but these are not critical enough to destroy the overall correlation. The results for water 
suggest that Rosenfeld-scaling will be useful for many other molecular fluids as a predictor of 
dynamical properties. 

We examine the diffusivity scaling with two different entropy estimators based on pair 
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correlation distribution functions (PDFs) . The first one which we refer to as the pair correlation 
entropy, S2, includes information from all atom-atom PDFs. The second one, referred to as the 
tetrahedral correlation entropy (Saa), considers only the contribution from the pair correlations 
between the tetrahedral atom-atom PDFs and is applicable to only water and ionic melts. We 
show that the isochore dependence of the diffusivity scaling parameters is most pronounced for 
Saa, much less pronounced for 5*2 and virtually negligible for Se- Thus the state-dependence of 
Rosenfeld scaling parameters increases as the structural approximants to the entropy become 
poorer. We demonstrate that isochores with similar temperature-dependence of the entropy 
estimators show very similar Rosenfeld-scaling parameters. 

The clear correspondence between the isochoric Rosenfeld-scaling parameters and the tem- 
perature dependence of the entropy along isochores has an interesting implication. The isochoric 
temperature dependence of the excess entropy of many atomic systems obeys a T~^/^ scaling, 
possibly because of the critical role of strong, short-range repulsions ^2 shows a similar 

7""2/5 dependence in simple liquids but not in anomalous fluids, as shown in this study and 
elsewhere |lipi9l H9]. This suggests that diffusivities for many atomic fluids will obey Rosenfeld- 
scaling best when the the thermodynamic excess entropy is used. The 2SRP fluid presents an 
interesting case in that Rosenfeld-scaling is applicable though there are deviations from T~'^^^ 
scaling. It would be interesting to examine these ideas in the context of recent work on strongly 
correlating fluids [5HI52]. 

Our results on structure-entropy-diffusivity relationships are experimentally testable. Atom- 
atom PDFs, diffusivities and calorimetric entropy are experimentally accessible as a function of 
temperature and density (or pressure). A comparison of diffusivity scaling with the thermody- 
namic entropy as well as the 0-0 pair correlation entropy between water and methanol should 
itself be very interesting, because they are both molecular, hydrogen-bonded liquids. Methanol, 
unlike water, does not show pronounced liquid state anomalies and we expect that this will show 
up in the diffusivity scaling with respect to the oxygen-oxygen pair correlation entropy but not 
with respect to the thermodynamic excess entropy. Since the atom-atom RDFs can be measured 
for essentially all molecular and polymeric fluids, our results suggest that excess entropy scal- 
ing should have substantial predictive value when combined with structural data, provided the 
complete set of atom-atom RDFs is used for estimating the entropy. However, coarse-graining 
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procedures which rely on mapping a subset of atom-atom PDFs will not necessarily preserve the 
thermodynamic state-point dependence of transport properties, specially in anomalous fluids, 
as has been noted in recent studies of coarse-grained potentials for water |46il47j . 

Our results suggest that Rosenfeld-scaling is valid for a wide range of simple and anomalous 
liquids. The deviations are small if the complete thermodynamic excess entropy is used, possibly 
because of the dominance of strong, short-range repulsions. A microscopic understanding of 
Rosenfeld-scaling is, however, necessary in order to understand the occurrence and magnitude 
of deviations in different systems. One possible route to developing such a microscopic approach 
is to use energy landscape approaches that have proved very useful in the context of supercooled 
liquids [531 154). Efforts to apply energy landscape analysis to understand Rosenfeld-scaling in 
simple and anomalous liquids have been relatively few |46 p55l[56] . 

Water-like liquids are systems where an additional length scale comes into play as a function 
of density and therefore deviations from Rosenfeld-scaling may be expected. As discussed in 
this paper, deviations are, in fact, small in this case. An analogous set of systems could be those 
where an additional length scale emerges as a function of temperature, as in inverse melting 
systems [57]. Some of the recent work on Rosenfeld-scaling in systems with finite repulsions, 
such as the Hertzian sphere and the Gaussian-core [58H61j . can be understood in terms of 
deviations from Rosenfeld-scaling behaviour due to change in the character of binary collisions 
as well as the possibility of additional length scales as a function of density. 

To summarize, the anomalous behaviour of the excess entropy (Se) and the associated scal- 
ing relationship with diffusivity is compared in five liquids with water-like thermodynamic and 
transport anomalies: SPC/E and TIP3P water models, Si02 and BeF2 ionic melts and the 
two-scale ramp (2SRP) fluid. By defining suitable length and energy scales for onset of anoma- 
lous behaviour, the density range of the excess entropy anomaly in water is shown to be much 
narrower in water than in ionic melts or the 2SRP fluid. While the Rosenfeld prediction of an 
exponential dependence of reduced diffusivity (D*) on the excess entropy (Se), D* = Aex.p{aSe) 
is found to be be almost quantitatively valid for all the systems, the exponential scaling pa- 
rameter, a, is found to show a small isochore-dependence in the case of water, presumably due 
to the rigid molecular structure. Unlike in the case of simple liquids, this isochore-dependence 
of the excess entropy scaling of the diffusivities is accentuated by replacing by the pair 
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correlation entropy, S2, for all the systems. In the case of the AB2 ionic melts and water, 
the effect is more pronounced if the diffusivity is scaled with respect to the pair correlation 
entropy associated with the tetrahedral atoms. More interestingly, variations in the the iso- 
choric temperature dependence of the entropic contribution (S'e, S2 and Saa) are found to be 
correlated with variations in the isochore diffusivity scaling parameters, suggesting that the 
organisation of the energy landscape changes significantly with density. The tetrahedral liquids 
(ionic melts and water) show very similar structure-entropy-diffusivity relations. Our results 
suggest that by systematically considering systems which represent different types and degrees 
of deviation from assumptions underlying Rosenfeld-scaling, interesting physical insights into 
the relationship between structure, thermodynamics and transport in liquids can be obtained. 
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Figures 

1. Plot of excess entropy with scaled number density, pj pm for (a) ionic melts and 2SRP 
fluid and (b) water models. The lines mark the respective minima and maxima in the 
2SRP fluid in part (a) and the SPC/E model in part (b). 

2. Correlation plots of Rosenfeld-scaled diffusivities with excess entropy, of (a) two-scale 
ramp potential (2SRP) fluid, (b) AB2 Ionic Melts (BeFa and Si02) and (c) SPC/E and 
(d) TIP3P. The lowest isotherms for BeF2 and Si02 are shown with filled symbols. Data 
points lying along the highest and lowest isochore for (c) and (d) are joined with smooth 
lines. The scaling parameter a for each plot is shown. In part (b), the datapoints at 
1500K(BeF2) and 4000K(SiO2) are shown as filled symbols, and are excluded from the 
fit. In part (c) and (d), data points having Se <-9.5 are excluded from the fit. 

3. Correlation plots of Rosenfeld scaled diffusivities, D* , in SPC/E water with: (a)S'2 (b)S'oo 
for selected isochores. Straight lines connect data from the highest and lowest density 
isochores. 

4. Different entropic contributions in SPC/E as a function of T~^/^: (a) thermodynamic 
excess entropy, Se, (b) pair correlation entropy, ^2 and (c) entropic contribution due to 

oxygen-oxygen pair correlations, x^Sqo with T"^/^. 
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SPC/E [30j 


TIP3P [31\ 


BeF2 133J 


Si02 [34J 


2SRP [16J 


p"^{g cm 


1.01 


0.98 


1.8 


2.3 


1.5 


r™(K) 


251 


195 


2310 


5000 


0.0548 



TABLE 1: State point corresponding to the maximum temperature along the TMD locus 
(p™, T™) for the systems studied. Potential parameters are given in the references. Values for 
the 2SRP are in the reduced units defined in the simulations details. 
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